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Radag^TrackiDR System 

This invention relates to radar tracking 
systems for missiles directed against airborne targets • 

The radar system may operate in an active mode 
in which the target is illuminated by radiation from the 
5 missile and in which reflected radiation from the target 
(termed "skin echo") is received by the missile during 
its flight. The reflected radiation is processed to 
obtain the bearings azimuth and elevation, and the 
speed or range of the target so that the missile can 
10 follow the changes of direction and speed of the target. 

An object of the present invention is the 
improvement of terminal accuracy in a pulse doppler 
radar s e eker . 

According to one aspect of the present 
15 invention, in a monopulse radar tracking system 

including either or both of a doppler tracking loop and 
an angle tracking loop, there are provided means for 
deriving from a target echo an intermediate frequency 
(I. P.) signal, means for estimating an I. P. target 
20 signal frequency periodically, digital filter means 
providing a frequency analysis of signals in the I. P. 
band, the characteristic of the filter „ means comprising 
a plurality of similar, ^eq^ntial, overlapping peaking 
characteristics, defirilng respective adjoining frequency 





bins, means for producing from target signals output in 
adjacent frequency bins^in the vicinity of the estimated 
target signal frequency a derived signal characteristic 
having a peak at a predetermined position in relation 
5 to the estimated target signal frequency, and means for 
correcting the estimate of target signal frequency and 
shifting the derived signal characteristic accordingly , 
the target signal resulting from said derived signal 
characteristic being employed in a said tracking loop* 
10 Where there is included a doppler tracking 

loop having a sf>eedgate filter in the I. F. signal path, 
the derived signal characteristic may have sections 
centred symmetrically above and below the estimated 
target signal frequency, means then being provided for 
15 comparing target signal constitutents of the respective 
sections and thereby determining the frequency error 
between the actual and estimated target signal frequency, 
the frequency error being employed to control the I.F. 
target frequency to tend to maintain it at a p>redeter- 
20 mined frequency within the pass band of the speedgate 
filter and to tend to bring the derived characteristic 
into alignment with the actual target signal frequency. 

In a missile-borne radar tracking system as 
aforesaid the width of the frequency bins is preferably 
25 controllable in dependence upon the estimated value of 
the difference between the target signal frequency and 
the predetermined frequency, within the speedgate pass 
band, the bin width being increased in response to a 
high estimated value to give good tracking ability and 
30 decreased in response to a low estimated value to give 
good velocity discrimination against targets moving at 
low relative velocities. 

In a radar tracking system including an angle 
tracking loop employing sum ajg# difference channels and 
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and means responsive^ to a ratio of the sum and difference 
channels and means responsive to a ratio of the sum and 
difference signals to* provide an indication of angular 
error between target line of sight and boresight, each 
5 of the sum and difference channels preferably employs 

means for producing a derived signal characteristic , each 
such derived signal .characteristic comprising a peaking 
characteristic similar to the individual characteristics 
of the digital filter means and centred on the estimated 

10 target signal frequency so as to provide a narrow , 
frequency pass band which tracks the estimated target 
signal frequency ♦ 

According to another aspect of the invention, in 
a monopulse radar tracking system employing sum and 

15 difference signals for the determination of target 

direction, the system comprising a doppler tracking loop 
maintaining an intermediate frequency target signal 
within a speedgate filter pass band, and digital filter 
means in each of the sum and difference channels providing 

20 at periodic update intervals ■ an analysis over a plurality 
of adjacent frequency bins of .target signal components 
within the said pass ban<^ and in a: method of confirming \ ^ 
the presence or absence of a target signal within a 
.-• particular frequency bin, a series of comparison processes 

25 are performed, each involving the sum of a predetermined 
number of successive power output values from the 
particular frequency bin accumulated with any previous 
such sums and a comparison of such cumulative sum with 
upper and lower threshold values which become progressi'v^ly 

30 closer with each comparison process, confirmation of the 
presence of a target signal in any comparison process 
being indicated when the uppers threshold is exceeded by 

the cumulative sum. confirmation oft. the absence of a 

**** 

35 target signal being indicated jwj&en the lower threshold 




exceeds the cumulative sum, and a further comparison 
process being initiated when the cumulative sum lies 
between the upper and lower thresholds. 

In such a method, a running average may be 
5 established incorporating a fixed number of bin power 
output values, the earliest incorporated bin power 
output value being discarded as a current value is 
incorporated so as to produce a running average bin 
power, this running average being compared with a 

10 predetermined threshold to provide a target-signal- 
present or absent indication the cumulative sum — 
indication and the running average indication 
contributing to a net condlusion in which the cumulative 
sum indication takes precedence if the running average 

15 indication points to a target signal absent, and in 
which the possible results of the cumulative sum 
indication: target signal absent, indeterminate, 
present, are treated as: indeterminate, present, present, 
respectively if the running average indication points 

20 to a target signal present ^ \ 
According to a further aspect of the invention, 
a radar tracking system employing sum and difference 
signals for the determination of target direction, 
includes digital filter means in each of the sum and 

25 difference channels providing at periodic update 
intervals an analysis over a plurality of adjacent 
frequency bins of potential target signal components, 
means for identifying a target frequency bin, means for 
applying sum and difference signals in respect of the 

30 identified target frequency bin as input signals to 
product means for producing a complex product of one 
of the input signals and the complex conjugate of the 
other, a signal-to-noise ratio indication being derived 
from the imaginary component of the complex product and 



a power level indication of the sum channel signal 
within the target frequency bin, the signal-to-noise 
ratio indication thus having a high value in the presence 
of incoherent reflections from multiple targets and a low 
5 value in the presence of coherent reflections from a 
single target, the system further including a basic 
signal-to-noise ratio indication derived from power level 
within the target frequency bin and average power level ^ 
over the plurality of frequency bins, this basic _ 

10 indication thus having a high value in the presence of . 
single or multiple targets within the target frequency 
bin and a low value in the presence of wideband noise, 
^ and means responsive to both of the signal-to-noise 
ratio indications to provide an indication of single j 

15 or multiple targets* 

.The relevant^ features of one embodiment of a 
pulse doppler radar seeker will now be described, by # - . 
way of example, with reference to the accompanying 
^ • . drawings, of which: : ^^^^ : ;^J^^ ^ -r - ^ '^*.\ x 

\ t diagram" form^ the basic /elements ~qf the '* r system , ^ V ^; i:^^ 

; ; ; Figure 2 showing a doppler tracking loop, and v ;1 

Figure 3 showing an angle tracking loop 1 
25 together with single/multiple targer discrimination | 
circuitry; "\ 



Figure 4 is a table illustrating a combination 
of target confirmation tests; 

and Figure 5 is a block diagram of an angle 
30 tracking and servo-control system. 
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In that part af the tracking system shown in 
Figure 1, an antenna 1 transmits radar pulses of 
frequency controlled by a local oscillator 4- and doppler 
tracking control path 6, the local oscillator and 
5 control signals being combined by a mix^r 8 and applied 
to the antenna by way of an amplifier 9 and circulator 
11, in known manner. The antenna includes a so-called 
comparator circuit which takes the target echo (and any 
other received signals) as received by a square array 
10 of four elements of the antenna and produces three 

* output channels, the azimuth and elevation difference - 

channels 13 and 15 and the sum channel 17, again in 
. known manner • - 

The antenna is steerable with respect to the 
15V missile by means of servo-controlled motors 19 to attempt 
to maintain the antenna boresight directed at a target* 
The actual direction of the boresight relative to the 
missile is determined by gimbal • pick-off • transducers 

20 ^ . r -^.-aa^e .sum and differeiwe signals are applied ^to t • 
\ - respective mixers 23, 25 and 2/ which each have an ; v 
; input from the local oscillator 4. /Che "transmitted ! ^1:, : 
frequency is the sum of the local oscillator frequency, 
a basic Intermediate (I. P.) frequency derived from a 
25 voltage controlled oscillator (VCO) 29, and the estimated 
doppler frequency ^corresponding to the voltage applied 
to the VCO 29 and being the output of the doppler- 
tracking loop. The output of the mixers 23, 25 and 27 
is therefore the IJ, plus the error in the estimation 
30 of the doppler frequency (true doppler frequency = f^). 

The mixers 23, 25. and 27*are followed by head 
amplifiers 33, 35 and 

Electronic angle tracing (E. A* T.) circuitry 

<A£fp6RFf 
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then follows in which a^pfespective proportion of the sum 

signal is added into each of the difference channels, 

the proportion being k for the azimuth difference and 

a 

k e for the elevation difference. This proportion is 

5 calculated to produce a net zero difference channel 

signal and thus simulate alignment of the boresight and 

target, the actual of f -boresight angle information (apart 

from error in the derivation) being contained in the 

control signals £ Q and £ (estimated bore sight /tar get 

a e 

10 . sightline error) applied to the E.A.T* , elements 39 and 

41 on lines 47 and 49. ; / , J \ - 2' : I \ 
The resulting sum and difference signals are " 
applied to range gates 53 » 55 and 57 which are open for ^ 
a predetermined short time at a controlled delay 
15 period after each radar pulse transmission. The timing 
of this delay period is controlled by a range tracking 
loop circuit 59* The sum channel is branched to feed an 
auxiliary range gate 51 which differs from the others in 
that the gated period is in two halves in one of which : . 
20 the gated signal is invert ed • ^ I .The output of , the jsumj^ I 
- I,, range gate 57 and the part inverted -output 'of ..the 

auxiliary range gate 51 are applied to ! a phase-sensitive-' 
. detector which thus gives a zero output if the received . 
signal lies equally in the first and second halves of the 
25 auxiliary gate. If the signal lies earlier or later, 
the PSD output will have a net negative or positive 
value which is used by the RTL circuit 59 to re-centre 
\ all of the range gates on the target signal. The target 
is thus effectively range tracked. 
30 The range tracking loop 59 also applies an 

inhibit signal to the transmitter amplifier 9 while the 
range gates are open to. ensure that no breakthrough of 
transmitted pulse through thej^r^ulator 11 can interfere 
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with target echo signals. ^ 

Following the^range gates 51, 53, 55 and 57, 
the sum and difference signals are applied to speedgates 
61, 63, 65 and 67 respectively. These are in fact band 
5 pass filters covering the possible range of doppler 
frequencies of interest and centred on the above 
mentioned intermediate frequency. The doppler tracking 
loop attempts to maintain the various sum and difference 
target frequencies in the centre of the speedgate pass 
10 bands by control of the VCO 29» ^ - ; - v ! ; / ■ ; 

- ■ - The sum and difference signals are then scaled 

"in accordance with the sum channel signal level by 
^ " age amplifiers 71 5 73, 75 and 77 controlled by an age * " 
detector 69* 

15 The signals until this point have all been 

analogue but are now converted by analogue digital 
converters 81, 83, 85 and 87 for the purpose of digital 
frequency analysis. ^In each of the sum and difference ^ 
channels is a f as t-f our ier- trans form (FPT) filter (91 f y : — 

20 93, 95 and 97 >^ 

;/ ; 1 of perhaps 32 , 128 , > 2^^f liter elements , : typically a - ; 

f \ ; say, 128, these filter elements respectively covering r 
;^ adjacent narrow frequency bands and together covering /v 
the speedgate passband. The frequency bands of such 

25 an FFT filter are known, and will be referred to in 

this specification, as frequency bins. The characteristic 
of each bin is of peaked form and overlaps with those of 
the adjacent bins. The bin outputs arise in digital 
form, representing the amplitude and phase of the signal 

30 component within the particular bin. One set of such 

output data requires information from 128 analogue samples. 
The output data rate of the FFT filters, i.e. the 
'update rate 1 is 78 Hertz i^# the FFT. overall bandwidth 
is 10 kHz and there are, 128 bins. Thus the update rate 

35 is equal to the individual b^^width. 



The output data rate of the FFT filters, i.e. the 
•update rate 1 may then be about 80 Hertz. 

The -width of the frequency bins is controllable, 
as indicated by control circuitry 89, thus covering the 
I. P. range with few, broad bins, or a greater number 
of correspondingly narrower bins. 

Each of the sunj, difference, and range channels 
is thus analysed into 128 (say) frequency bands within 
the IP doppler range and a target signal frequency can 
be identified very narrowly. 

Referring now to Figure 2, this shows the 
remainder of the doppler tracking loop, which is: entirely 
digital and is constituted by processes performed by a 
data processor. The input is derived on line 97 • from 
the sum channel FFT filter 97 in Figure 1, the single 
input 97* representing all of the 128 sum bin outputs. 
A target bin detection process 101 makes an initial 
estimate of the target bin identity, as will be explained, 



giving an initial target bin frequency f Q , (as a 
displacement from the speedgate centre frequency). 

A sliding filter process 105 to which the FFT 
filter-outputs are applied, produces, in response to the 
initial estimate of target frequency f c » a derived 
characteristic comprising two individual bin 
characteristics disposed symmetrically above and below 
f Q . One of these bins is derived from the left and 
centre ones of three adjacent bins and the other from 
the centre and right of the three. Discriminator (105), 
age (107) and divider (109) processes produce an error 
signal 6 d , being the frequency error between the estimated 
target frequency (f Q in this case) and the actual target 
frequency. 

The error signal £^jr is processed further, as 
will be explained, but in addition is directly added to 
the current target frequenc^ejs^timate f (initially 
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f Q ) in a process 111. The result is therefore the current 
magnitude of the : target doppler" frequency relative to the 
speedgate I.F. centre frequency. This digital value is 
subjected to -a gain control process G3 (113) and provided 
5 as a digital output of the doppler tracking loop. 

Referring again to Figure 1, the DTL output 
on line 113' is converted back to analogue form and 
applied to an integrator 121 which accumulates the doppler 
frequency error and controls the oscillator VCO (29) 

10 accordingly. With a constant frequency error t i.e. 

constant target acceleration, the loop will lock on to ; 

— the target frequency and cause the VCO frequence to ramp 
in tracking it. At constant relative target velocity 
the frequency error will cease, the VCO input will be 

15 zero and the WCO output will remain constant at the IF. 

An accelerometer 117 detects acceleration of . 
the missile and adds in a corresponding factor to the 
integrator input by way of adder circuit 119 • 

: Reverting to Figure 2, the further processing 

20 of the discriminator output, referred ,to above , is also " -7* 
employed to produce from a final addition 

f which is a further estimate of the target frequency 

... . - - — _ ' - ....... . ^ * ... . ' . . - * 

to speedgate centre frequency error. This signal f is 
employed in a target confirmation process 125 to be 
25 explained. 

Confirmation of the target bin is employed in 
several ways: the electronic angle tracking process 
referred to in Figure 1 is enabled; the guidance system 
of the missile is enabled; and incidental to the main 
30 arrangement a signal-to-noise indication is obtained 

for employment in multiple target detection, as will be 
explained. ' \ 

Considering Figure ^g&n greater detail now, 
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and in particular the target bin detection process 101, 
it will be recalled that the FFT circuits in each sum 
and difference channel provide outputs giving the complex 
amplitude in each of the FFT filter •bins 1 (typically 
5 16, 32, 64, 128 or 256 'bins 1 according to the bin width 
control setting 89). The filter bin width is equal to 
the total FFT bandwidth divided by the number of 
, filter bins and the output data rate is equal approx- 
imately to the bin width* 

10 Target detection is carried out on all of the 

power outputs (i.e. the square of the sum of the real ; * 
' and imaginary components) of the: sum channel FFT with : - 
the exception of the marginal bins at each end, i.e. " 
the first and last U L where N L - 9, 5, 3, 2, 2 for 

15 256, 128, 64, 32, 16 bin FFTs respectively. The power 
output of each bin is divided by the average of all the 
bins excluding the marginal first and last 1^, and then 
~ this value is compared to a threshold (T D ) which is , : 
determined so as to allow a certain number of crossings V 

20 if the input is pure thermal noise -("false alarms ) ^ S The 



technique used is J;o^ compare each ?bin ~* ^ 



the average power, as above, with the maximum scaled value 
of the previous bins and with the threshold" . V The J . • ■ 1 
maximum bin at the end of the process which also exceeds 

25 Tp is taken as a target alarm for the bin in question. 
If no bin exceeds the threshold then no alarm is found 
and the detection process is applied to the next set of 
FFT data. This detection method can be modified in 
practice to allow multiple alarms. 

30 If the target alarm is found* its bin number 

and centre frequency f is output so that the doppler 
tracking loop and confirmation processes can act on the 
correct frequencies. An additional calculation is 



performed in the detection process. The contents of. 
the alarm bin are reduced by the above threshold value Tjj. 
This is done because the probability distribution of . 
the maximum bin which exceeds a given threshold is 
5 approximately a constant (T^) plus a Eayleigh 

distribution (provided T D ~is > log e N). By subtracting 
the threshold value, the initial alarm can be treated 
in the same way as subsequent power outputs from the bin 
and this simplifies the confirmation system. 
10 The calculations performed in the above process . 

for each bin power p^ in order to establish the maximum ^ 
power bin are: ' : : ; v> * "r" "v 




15 where r^ is the scaled power of a bin numbered i; H is 

the total number of bins ; -is - the number of marginal i 

bins; and 3 is the h^her of a bin within those considered 

for averaging. If r ± T D and is also > r c max then ; 

";r c is set = r ± and i is stored as the alarm bin 

20 niamber. r^ max is the maximum value of r. 'so far 1 . The 
c 1 

stored value of i which remains when all bins have been 
considered is then the target alarm bin. This initial;, 
setting, or estimate, of the target frequency (f Q in 
Figure 2) is then employed to close the doppler tracking 
25 loop and initiate the target confirmation process. It 
will be appreciated that the target frequency, while 
being within the detected |>in (i), will not in general be 
precisely the bin centre^ frequency f Q . 

In order to get a feigh probability of target 
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acquisition with a low^ probability of 'false confirm* 
it is necessary to sum ^"Sbcoherently over a large number 
of FFT update periods* To make this possible, target 
tracking is necessary, since either the target must be 
5 kept in a fixed bin or the alarm bin must be known in 
some other way. For this reason a two stage acquisition 
process is used. -The first stage which has already been 
described detects the whereabouts of a likely target, 
the second, confirmation, ensures that this detected 
10 signal is indeed a target with a high degree of 
certainty. . . 

The method of operation of the confirmation :~ 
process, referenced 125 in Figure 2, is as follows:- 

**** 

the 'alarm 1 bin f is designated by the doppler tracking 

15 algorithm of Figure 2, except immediately after the first 
1 detect' when the detection subroutine provides this 
information. As in the detection algorithm, the 
equivalent of Equation (1) above is used to provide a 
scaled target/alarm bin output. (If the alarm bin is 

20 out of the range Nj+l to the scaled bin value is 

set to zero). Then this scaled 'target bin 1 value on 
one FFT data processing is added to the same quantity on 
the successive update. This summing continues until a 
predetermined number of FFT power output values from 

25 the target bin have been added together. At this point 
the sum is compared to two thresholds, an upper and lower 
threshold. If the sum of the bins is greater than the 
upper threshold then a 'target present' flag is set, 
if the sum is less than the lower threshold, 'target 

30 absent* is confirmed, or in other words the suspected 
target is rejected. If the £um is in between then the 
flag remains in a 'don't 'know' stkte. In this case a 
further set of 'target' bin poswer values from successive 
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FFT updates are addedVtV the f frst set and the accumulated 
total compared to two new and closer thresholds, the upper 
one again denoting 'confirmed presence 1 and the lower 
one • confirmed absence 1 . 
5 After a predetermined number of such comparison 

stages (provided as an input to the program) the upper 
and lower thresholds are made to coincide so that a 
definite decision is forced although it may not be 
conclusive, as will be seen. 

10 This successive comparison test is called the 

•cumulative sum 1 indication. 

At this point a running average of the last N 
sets, each of N FFT updates begins. As the current bin 
power value is incorporated in the running total the 

15 earliest value so incorporated is discarded. The 

criterion for confirmation is that the running average 
should be above a threshold (normally taken to be the 
last comparison process converged- threshold) and! 
that in addition the individual comparison processes 

20 should not reject. This latter test is put in to ensure 
reasonably rapid response if the target should suddenly 
disappear for some reason. 

The result is a 'running average 1 indication. 
The 'cumulative sum' indication and the 'running 

25 average 1 indication are combined as illustrated in Figure 
4 to provide a net conclusion. It may be seen that the 
cumulative sum indication takes precedence if the running 
average indication is negative, while the cumulative sum 
indication is, in effect, endorsed by one level of 

50 indication certainty if the running average indication 
is positive. 

In addition a certain apaount of re-initialisation 
< y ' j re- 

takes place. If the overall confirmation goes off the 
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•confirm present 1 stalte the. .running average is stopped. 
Also if the cumulative 5 sum test rejects or confirms, the 
accumulation of sub-set stuns ceases and each subsequent 
sub-set is tested separately. ^ 
5 When target confirmation (125) is complete the 

'confirm present 1 output 126 is used to enable the angle 
tracking loop and the missile guidance loop. 

The doppler tracking system employed allows 
doppler tracking to take place both within the digital 

10 processor of Figure 2 and through the analogue VCO 29 
in Figure 1. The key function that allows tracking to 
take place within .the data processor itself is the 
sliding filter/discriminator process, 103, 105, 107, 109. 
To form a frequency discriminator the power 

15 output of two neighbouring bins can be subtracted. If 

the target signal sits symmetrically astride the junction 
of the two bins the power in each will be the same and 
the difference will be zero. If the target signal is 
off centre one way or the other the result will be 

20 positive or negative accordingly. If two; such bins, 
centred on a specified target frequency estimate, can 
be simulated, therefore, frequency shift from this 
central position can be detected and the target signal 
frequency tracked continuously. 

25 To vary the position of the discriminator 

continuously in this way a technique of (effectively) 
sliding the filters of the FFT filter 97 is adopted, 
which requires the simulation of each of the two 
simulated bins mentioned above. The following applies 

30 to the simulation of each of the two. Given the complex 
output of two existing adjacent FFT bins it is possible 
to construct by using these two^qiian titles alone a new 
FFT bin which will have as^ itts^peak any chosen value of 



(4 



frequency between ttte^mid-po^ut of the two original bins. 
In order to derive a sj0£le algorithm a rectangular 
window is assumed and phase factors between adjacent 
bins of n /N are ignored* The new FFT bin characteristic 
can be written approximately as: 

f(X ) = a si ; 71 * + p si (i-x^ 1 " 3c) (2) 

'where x is the fraction of a bin from the centre of the 
first FFT bin where the response is required, a is^the - 
amplitude contribution of the first bin, and ^ is the 
amplitude contribution of the second bin* It can be 
shown that if a, 0 are related by a parameter 6 as 
follows : 

- 6 2 fl + 7tcotn6C 1-6)3 



a = 



7icotn6(6(l-6)(l-26)) - 1 + 26 - 26 2 



0 » 1 - a 

then the point f(x = 6) is a maximum of the filter of 
Equation (2)* ; Hence a sliding filter can be implemented* 
by taking the frequency f , evaluating the bins nearest to 
this and then calculating the parameters a, 0 needed to 
give two filter bins each separated by one half a bin 
separation" from the estimated target frequency (£). 
The same technique is used on the filters regardless of 
the type of window that is being used. The equations 
for the sliding filters thus become :- 

W ta = Integral part of ((f + 0.5B)t g +.1.5) 

x * -0.5 + N. - N - ft 

ta p s 

= a(x) (i~gj&. the above function a 

..... ^al^ulated for predetermined 
set of x values; 





C b - <1-C a ) b(H ta+ l) - Cl b(N ta ) 

where is the number of the bin 'intersected 1 by the 

required frequency, f is the estimated ^oppler frequency 
according to the digital tracking loop, B is the overall 
FFT bandwidth, t is the update time(« No, of binsfB)/ 
N p is half the number of bins in the FFT and b(N) is^ the 
the N tJl FFT bin (complex)* Thus C a and C b are the 
sliding filter outputs of simulated adjacent bins centred 
on the frequency f, and are again complex. 

It is then necessary to construct the 
discriminator 105 of Figure 2, and this can be done by 
taking the frequency error d = / C a| 2 ~t C b| 2# 
However, the sliding filter technique, results in a 
scaling at the origin of this quantity which depends on 
filter position and hence a correction for this has to 
be applied. A quadratic correction term is used of the 
form 

s d - a i + a 2 * + a e ^ ; - y ; „ ; i \ L 

where x = min (x, 1-x) and a^, a j take various values 

for the different window functions being used. This also 
allows the discriminator to be automatically scaled for 
the appropriate window function. 

The discriminator output 'd 1 must also be 
automatically gain controlled and this is accomplished 
by a simple first order feedback system operating off a 
square law detector 107 • The detector operates off a 
new bin of the form 

age ^a °b ^> 
which gives a filter centred on the target frequency. 
In conjunction with the disc^^inator this is believed 




to give optimum performance,. The age is thus 
agc n « agc Q-1 #+ (1-k) / C agc/ 2 S C 

where S is another scale correcting factor for the 
c 

different gains found at different positions of the 

sliding filter and different windows. This is given by 

— —2 
= 1 + x b n + x b~ 
c 12 

where b^ and b 2 are constants depending on the type of 

window, k equals exp (-t Vt„) where t „ is the age 

o age age 

time constant. Th.e resultant discriminator output £^ 
(from divider 109) is 

£d s d S c * > Hz 

The digital tracking loop is closed entirely 
within the data processor (see Figure 2) and provides 
the estimate of target frequency f that is needed for 
the confirmation process. The particular implementation 
used is a type two loop* A block diagram of the system 
is shown in Figure 2» 

It can be seen that there are two outputs of the 
loop, one for the confirmation process (f) and one 
feedback (f) to the sliding filter/discriminator. The 
transfer functions of these two are different and are 
given by 

J _ zG-^l+Gg) - G 1 

7 z 2 + z(G 1 (l+G 2 )-2) + (1-G^ 

7 A x + G 1 (G 2 -l)z 



7 + z(G 1 (l+G 2 )-2) + (1-G 1 ) 

Where f is the target frequency relative to the centre 
of the IF sum channel. These^ransf er functions are 



implemented by taking^he output of the discriminator £ d 
and then multiplying by^the gain (Gl). 

The algorithms for the implementation of these 
two closed loop transfer functions are as follows. If 
5 the discriminator output is £ ^ then 

f l (n) = *i(n-D + G l £ d 
f 2 (n) = f 2 (n-l) + G^Cn) 
f (n) = f-^n) + f 2 (n) 

f (n) = f-^n) + f 2 (n-l) 

10 These algorithms are illustrated in Figure 2 as 

follows. The output £ d of the discriminator processes 
105, 107, 109 is subjected to a gain Gl (127) producing 
an input G 1 £ d to a summing process 129 a storage delay 
device 131 feeds back the output of the summing process 

15 129 to its input at the following update. If f^n) is 
the output of the summing process 129 on the nth update, 
f^(n) must therefore be equal to G^ £ d V f^n-1), the 
loop therefore constituting a digital integration. 
The output of the loop, f^n) is applied to a further 

20 storage delay device 133 whose output f^(n-l) is applied 
to a summing process 135* A second input to this summing 
process is derived from the output f^(n) of the above 
loop as subjected to a gain G2 (137) to give G^-^n) and 
a further integrating loop 139, 14-1 to give f 2 (n) equal 

25 to f 2 (n-l) + G 2 f 1 (n). A further storage/delay device 14-3 
then gives f 2 (n-l) the other input to summing process 135. 

The output of summing process 135 is thus 
f(n-l) = f x (n-l) + f 2 (n-l) i.e. the estimated target 
frequency based on the preceding update of the PFT data. 

30 This output is applied to tl^#sliding filter process 103 




for calculation of the n^xi error £ d and is also 
applied to summing process 111 together with the current 
- error £ d as previously described, to provide the doppler 
tracking loop output. 
5 The signal f 2 (n-l) is applied to a further 

summing process 123 hut in this case with the current 
.input to process 135, i.e. f^Cn). The summed output of 
process 123 is thus f^n) + f 2 (n-l), denoted f and 
applied as the estimated target frequency to confirmation 
10 process 125* as previously described. 

Referring now to Figure 3 the FFT outputs from 
the azimuth difference filters 95 are input on lines 
... represented by 95* and the sum FPT outputs similarly 
on lines 97'. These sum and difference outputs, each , 
15 from, typically, 128 bins, are applied to respective, 
sliding filter processes 145 and 147 similar to those 
described for the doppler loop of Figure 2. In each 
— case^the estimated target frequency f derived from the / - 
/V doppier loop is used to select the bin in which it falls ^7 
20 ; a new or V supposed 1 bin being derived from this bin and 
, y y an adjacent bin, the^ or 
peak frequency , aligned with lthe ? estimated targe 
- Changes in the target frequency as determined by the 
doppler tracking loop cause the value of f to change, 
25 the simulated bin remaining locked to it and thus sliding 
up and down the I.F. range with the doppler error 
frequency and providing a high degree of target velocity 
discrimination. 

The sliding filter outputs are in digital form 
30 and are therefore applied to the digital equivalent of 
a phase-sensitive-detector. Thus the difference signal 
in complex form has its complex conjugate formed by a 
process 149, The sum and complex conjugate of the 



difference are then multiplied together by a process 151 
which produces real and^ imaginary outputs. The real part 
is divided by an age detector (157) output to give the 
amplitude ratio for each of the two channels i.e. 
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where S % the mean square of the sum signal, is the output 
of the age detector 157* Only the azimuth difference 
channel is shown but the elevation difference channel is 
treated similarly. This is fed to a Kalman filter in 
10 Figure 5 to give an output £ . £ is added to an 

aberration correction factor and fed to the EAT gain 
element to close the loop. The Kalman filter in Figure 
5 produces an estimate E* of the boresight error , that is, 
the estimated value of the target angle off boresight. 
15 This error estimate is combined in a summing process 

156 with a boresight direction indication relative to ^ 
the missile. The latter is. provided basically by gimbal 7" 
pick-off transducers 21 ( shown in Figured 1 also^^ 
corrected for radome aberration (150) which itself is : x 
20 dependent upon; boresight direction relative to the 
radome; radome temperature (163); and transmitted 
frequency (161). The aberration corrector thus produces 
azimuth and elevation angles of the effective boresight 
rather than the physical boresight. 
25 The resulting signal k a output from summing 

process 156 is employed in Figure 1 to control the EAT 
element 4-1 and close the EAT loop. 
A When the EAT loop is in equilibrium the value 

6 gives a true measure of J£ie boresight error. 
30 Considering now t£<e imaginary output of the 



product 151 1 this is divided by the age signal in divider 
153 to give a form of signal-to-noise indication. This 
output of the divider 153 gives an indication of the 
power of the signal components that are either incoherent 
5 between the sum and difference channels or in phase 
quadrature. The •incoherent* power component is large 
when thermal noise or jamming in either the sum or the 
difference channels is large , and tends to zero at high 
signal to noise ratio. The incoherent power output is 

10 also large on extended targets or on targets flying in 
formation that are unresolved by the seeker . 

The signal-to-noise (S/N) signal output from 
153 is applied to a dividing process 159 the other 
input of which is derived from the S/N output of the 

15 target confirmation process 125 of Figures 2 and 3. 
The two S/N indications differ in that the basic 
indication, from the confirmation circuit is derived 
as a ratio of the power in the target bin to the average 
power over the FFT band. It does not therefore 

20 distinguish between a coherent single target in the 
target bin and incoherent multiple targets within the 
target bin. TheQPSD S/N indication is derived solely 
from the target bin and distinguishes only between 
coherent (single) targets on one hand and multiple 

25 targets and noise on the other , giving a low value for 
the coherent target and a high value for the incoherent 
targets. If, therefore, the qpSD S/N indication is 
divided by the basic S/N indication from the confirmation 
process the result will be a low value for a single 

30 coherent target in the target bin and a high value for 
multiple targets in the target bin. The result of the 
division process 159 thus give^an indication of single 
or multiple targets. ^ p 




The system as described is adaptable for the 



various missile/targetwSbnditions that may arise in an 
engagement, and in particular in the terminal stages 
where conditions may change rapidly. 
5 The parameters in the doppler loop are made 

adaptive to allow for changing target conditions so 
that optimum performance is maintained at all times. 
Referring to Figure 2: 

(i) Immediately after target detection, gains Gl, 
10 G2 and G3 are set to high values so that the loop 
settles quickly. The gains are then reduced until, 
after a predetermined length of time, the gains take 
constant lower values to reduce the noise in the loop. 
The gains are increased as the missile flight progresses 
15 according to the signal strength measured in the age 
amplifiers (connection from age not shown), or by one 
or other of the signal to noise assessment indications 
described above, and according to the FFT bin width 
selected (see below). 
20 (ii) During the missile boost phase, the doppler 

tracking loop experiences substantial transients and the 
FFT bin width is selected to be a large value. Gains 
Gl, G2 and G3 are also set to high values to give a 
rapid response. 
25 (iii) If a target at long range is to be acquired, 

the FFT bin width is set to be. very narrow, and Gl, G2 
and G3 are set to low values after loop pull in. This 
optimises the acquisition threshold. 

(iv) If the target is at short range the FFT bin- 
30 width is set to be wide and Gl, G2 and G3 are set to 
high values. This optimises acquisition speed and 
minimises tracking errors. A wide FFT binwidth gives 
a faster update rate and hence gives the optimum 
performance very close to Ix^^t when rapid response is 

T 





important* Target racpr is determined from the range 
tracking loop (Figure 

(v) Once the loop has pulled-in, the signal f is 
proportional to the target frequency rate, i.e. the 
5 target acceleration. If f is large the FFT binwidth is 
set to a high value to allow a rapid tracking ability. 
The loss in doppler discrimination against other targets 
at similar speeds is not important since a manoeuvring 
target is not able to keep in close proximity to other 

10 targets. When f is small, a small FFT binwidth is used 
to give maximum doppler discrimination against targets 
flying in formation. 

(vi) When multiple targets in close proximity are 
detected (e.g. by the output from divider 159 in Fig. 3) 

15 or the target acceleration is small , the FFT binwidth 
is set to a low value to give maximum doppler 
discrimination. 

The overall angle tracking system can be 
conveniently divided into three parts . Firstly the EAT 

20 receiver, secondly the angle tracking filter and thirdly 
the servo control and stabilisation system. ^ 

A block diagram of the complete angle tracking 
and servo control system is shown in Figure 5* The 
angle tracking filter estimate c of boresight error 

25 angle in each of the azimuth and elevation channels is 
used to add a controlled portion of the sum channel 
which corresponds to these angles according to the 
stored D/S slope at the origin, into the difference 
channels, as described above with reference to Figures 

30 1 and 3« These difference signals then pass through 
similar IF chains as f or 5the sum channel, are sampled 
and analogue-to-digital converted (81-87) • They are 
then fourier transf 6rmed^S^-97>) to give the spectrum 
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of signals on the azimuth, and elevation channels. The 
sum and difference channels are divided by age signals 
digitally and elevation and azimuth phase-sensitive- 
detector outputs are formed (155) by taking the real 
5 part of D*xS; these outputs are proportional to the 
angular error between the true target boresight error 
and the best estimate of the boresight error. These 
calculations are all performed at the same update rate 
as the doppler and range tracking loops (i.e. at the 

10 update rate of the FFT output). 

After the PSD's (151) it is necessary to 
perform an axis transformation in order to transform 
from receiver axes to line of sight axes. This is 
implemented digitally. These signals are then passed 

15 through an angle tracking (Kalman) filter (171) which 
produces best estimates of sight line rate using a 
variable gain algorithm. The sight line rate signals 
r are used in various ways; firstly they provide an 
output to the autopilot as the guidance command, for 

20 which purpose they are transformed into missile axes; 
they are also used as the input to the EAT integrator 
which feeds estimate, of angle back to the. EAT element 
and finally it is used as input to the line of sight ' 
observer (173) to drive the antenna mechanism. This 

25 latter is added to an inert ial reference unit (IRU) 

gyro output (175) in order to provide the total movement 
of the target as measured in line of sight axes. 

These angular rates are still calculated at 
the (low) update rate of the loop. These rates are 

30 tfcen integrated using a high update digital integrator 
177 and the output is transformed from line of sight 
axes into motor shaft angles using digital axes 
transformations 179 • The angle Remands are then used 
to control a position controlj^oop ,181 for the antenna. 
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The angle demand 183 ?b subtracted froni the potentiometer 
pick-off s 21 and the er£or is fed into a controller 185 
which feeds demands to the antenna mechanism, servo and 
motors 187* The antenna/reflector plate is thus position 
slaved to the integrated IRU gyro and the target sight 
line error. A high update rate controller is used 
The effect of this is to remove body motions on the 
receiver output* An additional line of sight feed back 
189 is used to the EAT element (39, 41). This is added 
to the other feedback path and ensures that, the best 
estimate of sight line error can be used to subtract 
from the actual sight line error for EAT. The combined 
EAT feedback path has to be angle transformed (191) 
back into receiver axes and a radome aberration 
correction (159) is added in, as a function of the 
gimbal angles , as described in relation to Figure 3» 

The Kalman filter parameters are varied as a 
function of . conditions during missile flight and in 
general are arranged to give high bandwidth when the 
estimated time to impact is small. v^Timey to impact is :, 
estimated by dividing the missiie-target range (from 
the range tracking loop) by the closing velocity (from 
the doppler tracking loop). The filter bandwidth is 
also made a function of signal to noise ratio, the band- 
width being low when the signal to noise ratio is low. 

It may be seen that a versatile tracking system 
has been described. The sliding filters permit 
discrimination between closely spaced targets the doppler 
resolution being twice as good as with fixed FFT filters. 

An improved signal/noise ratio can also be 
obtained by sliding the simulated bin peak on to the 
target. ^ > 

Discrimination of ^multiple targets, e.g. by 
tracking one target and monitoring another, can be 



achieved merely by duplicating the sliding filters. 



